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ABSTRACT 


Advanced Ramjet Munitions Technology (ARMT) is an 


ongoing project within the Defence Advanced Research 


Projects Agency (DARPA) to research ramjet munitions.The ee 


e 
ARMT eight inch projectile uses ramjet thrust for a boosted 
trajectory, but operates on a thrust drag balance concept to 
create a pseudo vacuum trajectory during powered flight.The e 
trajectory was analyzed using an IBM 370 computer 
Simulation for three and five degrees of freedom. Work was ee 
also done to adapt the Ballistics Research Laboratories six eo 
degrees of freedom program to the IBM system.Projectile “ 
aerodynamic and mass properties were obtained from the : 
Norden Systema Wind Tunnel Data. Dispersion from the oo 
vacuum trajectory due to wind prior to ramjet burnout proved we 
minor.Dispersion due to constent thrust errors under 5% was au 
within a 600 foot radius at terminal guidance over a range .* 
of 33 miles. 
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I. INTRODUCTION 
A. ARMT OPERATIONAL CONCEPT 

Conventional artillery uses unguided projectiles to 
destroy a variety of targets including armored vehicles. 
Ballistic information and professional experience are used 
to direct fire onto a target with little chance of a direct 
hit in the first few rounds. A forward observer provides 
feedback to the guns for correction until the target is 
destroyed with subsequent rounds.Generally this approach is 
limited to the 30 kilometer range and often point target 
locations are not known. 

The ARMT projectile,formerly called AIFS, changes the 
conventional concept by providing a fire and forget,boosted 
and terminally guided projectile with a highly predictable 
trajectory.Conventional heat seekers and target identifica- 
tion systems are short range and require a blunt nose which 
has high drag. Active radar is not practical. Long ranges of 
60 kilometers are achieved by using solid fuel ramjet 
propulsion with an axigaymmetric air inlet spike. An 
accurately predicted trajectory is achieved by ensuring 
thrust equals drag to produce zero axial acceleration and 
create a paeudo-vacuum trajectory during powered flight. At 
a@ preset time (near ideal condition ramjet burnout) the 


inlet sapike ia diacarded and canardsa are deployed for 


ae 
4 


terminal guidance control. A passive or reflected target 
energy system will be used for terminal guidance. Research 
is ongoing and sponsored by DARPA. Projectile 


characteristics are at Appendix A. 
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Figure 1-1 ARMT Trajectory Concept 


B. THRUST DRAG BALANCE CONCEPT 

Balanced thrust drag implies that the thrust equals the 
drag vectors during powered flight. The aim of the balance 
is to negate meteorological variations throughout the 
trajectory by reducing all force vectors to zero except 
gravity. The resulting pseudo-vacuua trajectory is highly 
predictable. The balance is accomplished by measuring axial 
acceleration and maintaining it at zero.Appendix B explores 
thia concept in detail. The potential errors include 
transients due to wind induced reactions and measurement 
errors in the technical application of the balance in 
flight. Further, any asymmetries such as aerodynamic 
misalignment or thrust will produce a lateral force which 
must be offset by a small roll rate(sixth degree of 
freedom). This allows longitudinal forces to be the 
significant forces. 

The concept of thrust drag balance is best demonstrated 
in Figures 1-2 and 1-3.In Figure 1-2 the centre of gravity 
of the projectile is flying at an angle of attack (alpha) 
with a velocity vector (Veg) relative to a fixed earth 
reference frame. The aerodynamic forces are lift(L) and 
drag(D). Assuming static atatic atability and neglecting 


transient dynamics, the angle of attack becomes zero, which 
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reduces lift to zero.If thrust equals drag then only the 

b gravitational force exists. ’ 

: : If a crosswind(Vw) is applied, the projectile cocks 

into the wind because of the aerodynamic force acting on the oe 
centre of pressure. Adjusting the thrust and again assuming | 

static stability the drag can be offset as shown in Figure 

1-3 by applying the correct thrust.For future reference the 


component of drag normal to the plane containing Vcg is Dw. 


ome 
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Lift as 


Drag 


Drag Thrust 


Figure 1-2 Projectile Dynamics Thrust Equals Drag 
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Figure 1-3 Projectile Dynamics With Crosswind 


If a tailwind or headwind ia applied, thrust is adjuated 
to compensate in the same manner. 

The dynamics or transient effects due to wind which 
create damped oscillatory motion have a cumulative net lift 
effect( as shown in Figure 1-4) aa angle of attack is 


proportional to lift.The area under the curve in Figure 1-4 


18 


does not average to zero due to the damping. For future 


reference thia lift due to wind will be referred to as Lw. 


Figure 1-4 Angle Of Attack Versus Time 


Calculations integrating alpha(t) twice to demonstrate the 


wee 


above in the yaw or pitch plane are at Appendix B. 


Cc. THESIS SCOPE 

It is emphasized that the scope of this thesis is to 
analyze dispersion from the pseudo-vacuum trajectory at 
terminal guidance time for the eight inch ARMT projectile.In 
colloquial terms the “basket “ from which the projectile 
will be guided to target was determined. Further, this 


basket is based on thrust errors and wind effects for a 


13 


launch angle of 45 degrees and muzzle velocity of MACH 2.2 
at sea level. Linearized aerodynamics were used in all 
calculations. No attempt was made to establish a “footprint” 


for projectile impact. 


II. TRAJECTORY ANALYSIS METHODOLOGY 

A. VACUUM TRAJECTORY ANALYSIS 

A short and simple BASIC program was written to 
establish the vacuum trajectory and fuel burn time of the 
projectile in a two dimensional scenario where thrust 
exactly equaled drag. The Standard U.S. Atmosphere was used 
and the drag coefficient linearized [(Reference 1]. The 
program listing and baseline data are at Appendix D. The key 
parameters and equations used are listed at Appendices A and 
D.The baseline pseudo-vacuum trajectory is shown in Figure 


2-1. 


B. THREE-DEGREES-OF-FREEDOM MODEL 

The computer simulation was established using the 
Interactive Ordinary Differential Equation (IODE) library 
program at NPGS. A listing of the three degrees of freedom 
equations is in the IODE specifications at Appendix E. The 
Simulation is based on wind tunnel results from Norden 
Systems (References 2 and 31] and assumes a horizontal plane 
at 30000 feet altitude.This height was chosen as significant 
flight time is at altitude and the time of flight in this 
plane was near actual flight time. Primarily the model was 
used to study crosswind effects but resonance with the 


projectile short period was simulated. Different thrust 
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scenarios were also used. Figure 2-2 depicts the model 


f the X-Y plane. 


- burnout time=101 secs 
i burnout range=175100 feet 
burnout height=10990 feet 
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Figure 2-1 ARMT Pseudo-vacuum Trajectory 


rmirrwe 


Y 
(inertial) 
mx =Lsiny — Ocosy 
my = Lcos ¥+ Dsin 


fy (body } 


X (inertial ) 


—-™ x(body) 


Figure 2-2 Three Degrees Freedom Model (Yaw Plane) 


The crosswinds were varied in intensity and duration to 
analyze crossrange dispersion. As well thrust was varied as 
a function of drag by scaling the drag. Cmq predictions were 
approximately -190 to -205 by Norden Systems (Reference 2] 
and an analysis was done to determine the effects of varying 
this damping coefficient. To further study the projectile 
characteristics, trajectories were simulated with no thrust 


or mass attenuation with varying wind profiles. Figure 2-3 , 


fit 
Gs 


16 + #(Cmg)6 + £(Cmg)@ = #(Cmg) ¥ 


shows the baseline three-degrees-of-freedom trajectory of S 


thrust equals drag with no wind. 
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Figure 2-3 Baseline Three~Degrees-Of-Freedom Trajectory 
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C. FIVE-DEGREES-OF-FREEDOM MODEL 

Again the IODE program was used to solve a five degrees 
of freedom system of equations. The roll plane was neglected 
for the following reasons: 

-the program could not accept sufficient functions to 
define six degrees of freedom. 

-the fins are expected to damp out initial spin due to 
gun rifling (900 RPM) in the first few seconds of flight. 
Sufficient spin to average projectile asymmetry errors to 
zero but not enough to significantly affect dynamics will 
remain.(ie. gyroscopic effects and Magnus effects) 

-the initial spin damping dynamics are beyond the scope 

of thia thesia. 
The projectile was assumed to be axisymmetric allowing for 
ease of calculations in the yaw and pitch planes. Figures 
2-2 and 2-4 depict the yaw and pitch planes for the 
model,respectively. 

A listing of the equations is provided in the IOQDE 
specifications at Appendix F. Trajectory simulations were 
run with varying wind profiles and thrust error scenarios. 
Terminal guidance time (Tc) was determined by the baseline 
five degrees of freedom trajectory shown at Figure 2-1 as 
101 seconds. This trajectory has no wind and thrust exactly 


equaled drag until discard of the air inlet spike and 


fis 
ca 


eanard deployment.All aimulationa aasumed a 45 degree launch or 
angle and Mach 2.2 muzzle velocity. Dispersion from the 
pseudo-vacuum trajectory was measured at terminal guidance oe 
time including the following variables: 2 

-X=range at 101 seconds 

-Y=crossrange at 101 seconds 

-Z=height at 101 seconds 

-Tb=burnout time 

-Vcg=speed at time 101 seconds 2 


Dispersion is defined as the difference between value for 


perturbed trajectories and the baseline trajectory. 
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Figure 2-4 Five-Degrees-Of-Freedom Pitch Plane Model 
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D. SIX-DEGREES-OF -FREEDOM MODEL 

The Ballistics Research Laboratories (Aberdeen) six 
degrees of freedom model HTRAJ was obtained and adapted for 
use on the NPGS IBM 370 system. The adapted listing is 
included as Appendix G. Program compilation was achieved but 
output data was not satisfactory due primarily to input data 
syntax problems. Further,the program requires extensive 
modification to be able to vary thruat asa a function of drag 
with the ease of the previous models already described. Time 
was not available to correct these problems. Detailed 


documentation on the program HTRAJ is available through BRL. 
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III. TRAJECTORY RESULTS 


eee ee See OS SS eS 


A. THREE DEGREES OF FREEDOM 
1. Crosswind And Thrust Error Sensitivity aS 

Figures 3-1 and 3-2 show the effects of constant 
thrust errors as a percentage of drag with no wind.Recall 
that the three-degrees-of freedom trajectory is at constant 
altitude.The figures therefore represent a thrust error 
sensitivity in terms of burn time and range variations due 
to drag or thrust.The trajectories all pass nearly through 
the 218000 foot range mark at t=100 seconds, but the burnout 
ranges are significantly different due to the varying burn 
times resulting from projectile thrust.Figure 3-2 
demonstrates the obvious conclusion that the burnout time 
decreases as thrust increases. 

Figure 3-3 combines the effects of constant 
crosswinds and thrust errors to demonstrate resultant cross 
range dispersion.The constant wind velocity lines are very 
nearly linear, and the pattern is essentially a shift about 
the zero-zero point. If thruat is leas than drag, the 
diaperaion ita in the direction of the wind vector.If thrust 


ia greater than drag then the overcompensation producea a 


dispersion opposite to the wind vector. 
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Figure 3-1 Thrust Error Sensitivity (Range) 
Three Degrees Of Freedon 
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Figure 3-2 Thrust Error Sensitivity (Burntime) 
Three Degrees Of Freedom 
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Figure 3-3 Crosawind Sensitivity Three Degrees Of 
Freedon 
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Of interest was the dispersion due to the same wind - 
profiles without thrust or mass loss. This represents an 
absolute dispersion with a given wind profile. These results 


| are in Table 3.1. 


a Table 3.1 No Thrust Crosswind Effecta 
. Three Degrees of Freedom 


Thrust=0O t=125 seconds 


t 

' vw (ft/sec) dispersion (feet) eas 

crossrange 7s 

a fe) fe) 

- 60 195.76 7 
90 293.80 

ue 120 392.02 

o 

minimal effect in the thrusted scenarios.The dispersion in 


: The data clearly shows that crosswinds have a 
the thrust-equals- drag case is due to the transients (Lw) 
> 


and is negligible as well. This is illustrated at Table 3.2. 
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Table 3.2 Dispersion Due to Net Lift 
from Transients 
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Thrust=drag 


vw (feet/sec) dispersion (feet) 
crossarange 
te) re) 
60 -018 
3930 -027 
120 -036 
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Angle of attack damped to virtually zero during the 
trajectories (see Figure E-3) but the absolute maximun 
values were also recorded.Figure 3-4 demonatrates the 
relation between angle of attack (alpha) and the wind 
velocity.Thruat errors had virtually no effect on the angle 


of attack magnitude. 
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Figure 3-4 Wind Effects On Maximum Angle Of Attack 


2. Damping Coefficient Effects 
The magnitude of Cmq was determined to be 
approximately -200 in the flight speed ranges of the ARMT 
projectile {Reference 4]. To determine if this coefficient 
Significantly affected dispersion Cmq was varied from 0 to 
1.9 x 106 and the dispersion recorded. Figure 3-5 shows the 


relation between Cmq and dispersion (crossrange) with a 


constant wind of 60 feet/second and time of 125 seconda. 
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Figure 3-5 Damping Coefficient Effects 
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Careful examination of Figure 3-5 shows that dispe- 
rsion increases slightly for realistic values of Cmq and is 
therefore proportional to Cmq.Due to the small sensitivity 
of dispersion to Cmq for all practical purposes one can 
conclude dispersion is independent to Cmq.(The increased 
disperaion proportional to Cmq can be explained as the 
difference in the areas under a highly damped sinusoid 


versus a non-damped sinusoid is considerable.) 


3. Sinusoidal Wind Analysis 

To test the effects of phase angles and wave- 
lengths on the projectile dynamics a sinusoidal wind was 
input to test the model.The wind (Vw) was of the form Vo 
sin(2 PI Vo t/lambda + phase) where Vo was set at 60 
feet/second.Figures 3-6 and 3-7 show the trajectories for 
the sine and cosine winds to burnout at thrust equal to 90% 
drag.(Lambda=range and lambda=.5 range)Of note is the fact 
that the sine wind produces a continuous dispersion while 
the cosine wind does not when thrust does not equal 
drag.This discovery led to an expanded study of different 
Phase and wavelength sinusoidal wind inputs for various 
thrust acenarios.As well a simplified model was constructed 
as shown at Figure 3-8 and calculations confirmed the three 
degrees of freedom model results.Figures 3-9 and 3-10 depict 


the results of the calculations. 
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Figure 3-6 Sinusoidal Wind Effects 


i) 
NJ 


110 120 


(Sine) 


130 


140 


ico] 
a 
nN 


1.38 


1.16 


DISPERSION (FEET) 
0.82 


OF a taie aie ee aS aie OS 6 Saas 618 SS iS SS ee ea wo a 6 Sa a Se oe a ee ee oe 


0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 
TIME (SECONDS) 


Figure 3-7 Sinusoidal Wind Effects (Cosine) 
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Figure 3-8 Simplified Three Degrees Of Freedom Model 


The simplified model was used to analytically 
evaluate the crossrange dispersion (y) in the cases of asaine 
and cosine winds to confirm the validity of the three- 
degrees-of-freedom model.Calculations follow: 

General equations of motion in the y direction fron 


Figure 3-8 are: 
M dVy/dt =Dy which approximately equals D(Vw-Vy)/Vxo 
if Vxo=Vo=constant.Then transforming the equation we get: 


dVy/dt + D/(MVo)Vy = D/C(MVoO) Vwlt) 
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or in time differentiated notation: 


Y + By = BVw(t) if D/CMVo) =B 


If a cosine wind Vw(t)=Vwo cos (2 PI Vo t/A) is applied 


and the differential equation solved we have: 


y (t/a =K le-Bt -cos(2 PI Vo t/r) + 


DAsin(2 PI Vo t/A) 
2 PI m Vo2 


let K= 2 1/2 m Vo2 Vwo 


Figure 3-9 graphically depicts the results as t increases.A 
zero asymptote is observed and this conforms with model Bye 
results. 


Now if a sine wind Vw(t)= Vwo sin (2 Pi Vo t/A) is applied 


and the differential equation solved: aes 
y(to/A =K DX €1-cos(2 PI Vo t/~)) - 


4 PI 1/2 m Vo2 
sin(2 PI VO t/a) +4 _PI(1/2 mVo2) (1-e7Bt) - 
DA 
Figure 3~10 graphically depicts these results as t 
increases.The asymptote is 


an K | (4 PI 1/2 m Vo2) (1-e7-Bt) 
“. D> 


and this conforms to the model also. 
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Figure 3-9 Cosine Wind Analysis Results 
41 
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Figure 3-10 Sine Wind Analysis Pesults 


The constant wind case was also derived: 


Vwlt) Vw and the solved differential 
equation is: 


y(t)= Vw t - Vwm Vo/D e-Bt 


and as t increases y approachea Vw t. 


Specific calculations were done to confirm 
dispersion results at burnout.The calculations were based 


on: 


y=dy/dx X = 1/Vo dy/dt X 

Vwo=60 ft/sec 

Vo=2190 ft/sec 

A =304000 ft 

t= ~/Vo 

D=1/2 9 Vo2 A Cd (=226 lbf) 
Thrust=.9D therefore D=22.6 lbf 


M=225 Ibm 


Taking the time derivative of the cosine wind y(t) 


dispersion formula derived earlier we get 
dy/dt=K|-Be-Bt + 2 PI Vo sin(2 PI vo t/a + os 
A en 


DA , 2 PI Vo cos(2 PI Vo t/a) 
4 PI 1/2 mVoz » 


ow aa slo, xh 
7 rian ieee era 
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at t= A/Vo the equation becomes: 
dy/dt = K(-Be-BA/Vo + B) 


and substituting the values into the approximationn 
for ys: y=0 ft. 
Now, performing the same operations on the sine 


wind dispersion equation : 


dy/dt=K DA . 2.PI Vo sin(2 PI Vo t/”) 
4 PI 1/2 mVoz ~ 


-2 PI Vo cos(2 PI Vo t) +(4 PI 1/2 m Vo2), 
x x DA 


(Be~Bt) 


and at t=WVo 


dy/dt =K] -2_PI Vo -4 PI 1/2 m Vo2-Be-Bt 
» DA 
and y=20 £t approximately. 
Finally, taking the constant wind dispersion 


formula and taking its derivative ; 
dy/dt = Vw- Vw e~Bt 


and y= 10 ft approximately. 
These results are consistent with the three-degrees-of- 
freedom model conaidering the initial approximations. 
To study the short period resonance phenomenom the 


Phase and wavelength wind variations were documented in 
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detail for a condition of thrust equal to 90% drag.Figures 
3-11 to 3-13 show the crossrange dispersion effects from 
wavelengths of SOOO ( lambda = .016 range) to 
303000(lambdatrange) feet and the phase angle variations 
from OO to 2 PI.At approximately 9850 feet the dispersion 
characteristics were no longer symmetric but remained under 
one foot.This is demonstrated at Figure 3-14.A crossrange 
dispersion contour was constructed at Figure 3-15 but this 
did not emphasize the resonance phenomenon sufficiently due 
to linearized aerodynamics and small dispersion 
magnitudes.The parameters that varied most significantly 
near resonance were angle of attack and burn time.Indeed,it 
is clear that the projectile is unstable in this region with 
a maximum angle of attack of 1.21 radians(70 degrees).Due to 
the increased angle of attack, drag increased and burn time 
decreased in the order of seven seconds (5%) compared to the 
baseline trajectory. This clearly affected the range at 
burnout .Due to large angle of attack the linearized 
serodynamics are not valid.Projectile performance at wind 
resonance is worthy of additional analysis.It was decided to 
demonstrate resonance by plotting the maximum angle of 
attack (alpha) versus wavelength with different thrust/ drag 
scenarios.The resonance results are at Figure 3-16.The 
projectile natural frequency is estimated at 0.22 Hz. 


Resonance was confirmed in the 0.22 Hz regime for thrust 


Bi 


errors under 10% including varying Vwo.This was necessary 45 
the system is non-linear and Vwo could have affected 
resonance.It was noted that the curve shifted and increased 
or decreased in magnitude slightly in each case.This is due 
to the non-linearity of the system equations and the 
variations in Vo during flight.The absolute c.se of no 


thrust showed resonance at 10000 feet with a maximum angle 


‘of attack of 1.25 radians(72 degrees). 
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Figure 3-11 Dispersion Versus Wind Phase Angle(Large Lambda) 
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Figure 3-12 Dispersion Versus Phase Angle 
(Medium Lambda) 
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Figure 3-14 Dispersion Versus Phase Angle 
Resonance Lambda 
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Figure 3-15 Phase Angle Versus Wavelength 
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Figure 3-16 Detrimental Resonance Between Projectile 
Short Period and Spatial Frequency of 
Wind Profile 
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It is clear that if the wind profile experienced by 
the projectile has harmonic content near the natural 
frequency then the round will be unstable and/or have large 
dispersion at burnout. This is best illustrated by viewing 
the wind as Vw(t) with reference to the projectile. By a 
Fourier transform of the wind profile resonance with 
projectile natural frequency (fn) will be evident in the 
frequency domain.The wind oscillations need only be velocity 
variations with the appropriate harmonic content.Figure 3- 
17, which illustrates this concept, shows fn coinciding with 
@ maximum of Vw(f). One would anticipate stability problems 


for this case. 
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Figure 3-17 Fourier Transform Of Wind Profile Concept 
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Calculations for the long period resonance or 
phugoid show that the phugoid is irrelevant as flight time 
is too short: 

T=4.44*Vo/g=301.8 seconds (Reference 3] 

The results of the sinusoidal winds were also 
compared with the aerodynamic jump phenomenom outlined in 
BRL report 1077 (Reference 4 1. At first glance the 
Sinusoidal wind results appear consistent with aerojump,but 
detailed calculations show otherwise. The rate of change of 
initial angle of attack is critical to aerojump and is 
easily incorporated into the aerojump equation.The results 
for the sine wind are a function of wavelength and of the 
same form as the simplified model results but the cosine 
wind and constant wind both yielded a small aerojump 
proportional to the initial wind.The magnitudes of the 
aerojump dispersions were more consistent with the net lift 
transient effects. Further, continuous perturbations in the 
form of wind after time zero are not in the aerojump 
equations and for all practical purposes the phenomenom is 
not observed in the model.These results should have been 
apparent before calculations were performed. Aerojump 
equations are not appropriate for non-zero wind 


conditions.Clearly, when wind is included the aerojump 


equations cannot predict dispersion. 
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B. FIVE DEGREES OF FREEDON 

As previously stated dispersion was measured from the 
ideal pseudo-vacuum trajectory at a fixed time for 
initiation of terminal guidance (Tc) 101 #£seconds.Thrust 
error sensitivity in the zero wind condition is presented in 
Figures 3-18 and 3-19.In these two figures it is clear that 
thrust errors can have a significant impact on the location 
of the projectile at the Tc.Figure 3-18 indicates that as 
thrust increases the dispersion moves downrange as 
expected.Figure 3-19 shows that as thrust increases’ the 
height of the projectile at Tc decreases. Due to added 
thrust the projectile is located at a point further along 
the parabolic trajectory; hence the height is less. 

As well Figure 3-20 shows these effects on available 
burn time of the ramjet.For thrust equals drag burnout at 
101 seconds exactly equals the Tc.With less thrust, terminal 
guidance is initiated before all the propellant is burned. 
Conversely, when thruat exceeds drag, burnout occura prior 
to initiation of terminal guidance.Assuming this fixed time 
for initiation of terminal guidance, the projectile may 
experience up to seven seconds of unpowered flight ¢ 7% 
error) before the guidance is deployed. The period of 
unpowered flight with associated high drag will cause major 


dispersions as shown in subsequent figures. 
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Figure 3-19 Thrust Error Effects (Height) 
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Figure 3-20 Thrust Error Effects (Burn Time) 
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Crosswind sensitivity results were very much the same 
as the three degree simulation.In fact the results. appear 
even more linear with the same rotation about the zero zero 
point as in Figure 3-21.The gross features of Figure 3-3 and 
3-21 are identical but Figure 3-3 ashowa diapersion at 
burnout whereas Figure 3-21 shows dispersion at Te. 

The crossrange dispersion results of a no thrust 
projectile with constant crosawind profiles and with thrust 
equal to drag (transients effects) are at Tables 3.3 and 3.4 


respectively. 


Table 3.3 No Thrust Wind Effects 
Five Degrees 


Thrust =0 time=101 seconds 


Vwxy (feet/sec) diaperaion (feet) 
crossrange 
fe) 0 
30 92.51 
60 185.11 
90 272.90 
120 370.97 
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Figure 3-21 Crosswind Effects Five Degrees 
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Table 3.4 Dispersion Due to Net Lift 
from Transients 


ee eee 


Thrust=drag time=101 secs 


Vwxy (feet/sec) dispersion(feet) 
fa) fe) 
30 ~0202 
60 -0413 
30 -0643 
120 -0919 


Burn time was not significantly affected by crosswinds. 
Headwind and tailwind effects on range,height and 
burntime in different thrust scenarios are shown in Figures 
3-22 to 3-26. It is clear that dispersion from the pseudo- 
vacuum ivayecteey is minimized when thrust nearly equals 
drag.It is also apparent that the overthrusted trajectories 
involve greater dispersion as burnout occurs prior to 
terminal guidance and unpowered flight occurs for some time 
dependent on the scenario.This same error is observed in the 
case of headwinds and is most clear in the Figures 3-22B to 
3-25B about the thrust equal to drag axis. (Note that Figure 
3-22 to 3-25 inclusive have A,B and C sub-figures which are 


simply enlargements of key portions of the respective main 


figure). 
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Figure 3-22 Headwind Sensitivity (Range) 
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Figure 3-22A Headwind Sensitivity (Range) 
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Figure 3-22B Headwind Sensitivity (Range) 
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Figure 3-22C Headwind Sensitivity (Range) 
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Figure 3-23 Tailwind Sensitivity 
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Figure 3-23A Tailwind Sensitivity 


ee 


LEGEND 


(Range) 


Veit, elas sia Tere ata 


Pi ae ee ae 


Stes ey De 


OSESLI 


ooesal OSescZi ooesal OSTISAI OOTSLI ososcZi OoOoscZI OS6FLI eee Oks 
(Ladd) ADNVY 


=) 
Z, 
fx] 
Oo 
fx) 
4 


oF) 

2 

° 
io) 
- 

nm % 

356A 
gE 
aw 

oO 

om 
es 
ica) 

on 

S35 
me 
mm 
e 

= 

r 

iis 

ri 

° 

o 


' 


~~ 


Figure 3-23B Tailwind Sensitivity (Range) 
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Figure 3-23C Tailwind Sensitivity (Range) 
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Figure 3-24 Headwind Sensitivity (Height) 
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Figure 3-24A Headwind Sensitivity (Height) 
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Figure 3-24B Headwind Sensitivity 
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Figure 3-24C Headwind Sensitivity (Height) 
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Figure 3-25 Tailwind Sensitivity (Height) 
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Figure 3-25A Tailwind Sensitivity (Height) 
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Figure 3-25B Tailwind Sensitivity (Height) 
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Figure 3-25C Tailwind Sensitivity (Height) 
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Figure 3-26 Tail/Headwind Effects (Burntime) 
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Deviations from the ideal trajectory were correlated 
and tabulated in Table 3.5.Results were conservatively 
linearized and are presented in that manner. Wind errors are 
read in addition to the thrust errors at the left of the 
table.For example: 


If thrust error=-5% 
crosswind=120 ft/sec 


predicted dispersion is X=-540 ft 


Y=19 £t 
2=130 ft :from ideal basket. 


Table 3.5 Dispersion Analysis Results 


Thrust Errors Winds 
continuous 120 ft/sec 


Cross Head Tail 
Errors Y xX Zz 


Oo %* 


f Velocities at time 101 seconds were measured and 


are 


shown at Figure 3-27.The total variation of extreme cases is 


from 2285 feet/sec to 2315 feet/sec, 
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Figure 3-27 Velocity Sensitivity 


or only 30 feet/sec. 


Wind effects during near thrust equals drag portions of 
the trajectory were very small.The critical factor 
affecting the dispersion will be the setting of the terminal 
guidance clock in the projectile.This presents several 
options for setting initiation of terminal guidance: 

~Fixed time irrespective of trajectory Cie. 
elevation of gun).Burn time would be shorter at less than 45 
degree launch and the basket would be more sensitive to wind 
and thrust errors.As well the trajyectory would change and so 
would the predicted basket. 

-Fixed time for a given trajectory.This requires 
analysis of the sensitivity to elevation of the gun and a 
manual Tc setting method for gunners. 

-Ramjet burnout.This would ensure that Te occurred 
at an “optimum” point almost exactly on the pseudo-vacuum 
trajectory but the location along it at Tc would be 
difficult to predict. 

-At given pressure for altitude (Reference 5 ].This 
would establish accurate baskets for given gun elevations 
assuming the altitude chosen was before ramjet burnout. 

Several trajectories were simulated with a varying wind 
profile as a function of altitude in the Standard 
Atmosphere.These are shown at Figure 3-28 as compared to the 


pseudo-vacuum trajectory. 


8 


Figure 3-28 Three Dimensional Trajectory 
Dependent Wind Profile 


IV. CONCLUSIONS 

Generally this thesis confirms the viability of th 
pseudo-vacuum trajectory from ae thrust equals dra 
projectile.Indeed winds up to 120 feet per second and thrus 
errors up to five percent created dispersions of less tha 
600 feet over a range of 175090 feet for the ARM 
projectile. 

Specifically the following conclusions are made: 

-the drag due to wind (Dw) is the aerodynami 
force which creates the most dispersion. 

-the net lift due to transients (Lw) is a mino 
serodynamic force in creating dispersion. 

-the dispersion is proportional to the dampin 
coefficient (Cmq),although for practical values of Cm 
dispersion can be treated as insensitive to Cmq. 

-the projectile resonance or short period is a 
approximately 0.22 H2 and is unstable in th 
resonance regime due to variable wind profiles.The phugoi 
is irrelevant. 

-thrust errors are the most significant error 
contributing to dispersion. 

-windsa have a minor effect on dispersion durin 


powered flight when thrust approximates drag 


except where harmonic Content near resonance is 
contained in the wind profile. 

-headwinds create greater dispersion as ramjet 
burnout occurs before terminal guidance is 
initiated.Crosswinds and tailwinds generally 


have minimal effects. 


Recommendations for further study related to 
thesis are: 
-Dynamics at launch including the initial 
damping, Magnus effects and aerodynamic jump. 
-Dynamics at terminal guidance including air inlet 


spike ejection and canard deployment. 


-Complete adaptation of BRL HTRAJ on the NPGS IBM 


system to use non-linear aerodynamics. 
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Figure A-1 ARMT Projectile Baseline Design 


The following characteristics are obtained from 
References 2 and 6. — 
length=72 inches : 
diameter=8 inches 
launch weight=225 lbm - 
launch CG=39.4 inches from spike apex pis 


fuel weight=27 lbm (CG at 52.0 inches from apex) oat 


84 ae 


impul 


Mass 


Aerod 


s8e=950 “seconds” 

Property Equationa: 

CG=(225x39.4)-(52.0xfuel burned)/(225-fuel burned) 
Izz=.35-¢(.0019xfuel burned) in slug-ft 


Ixx=TIyy=18.8-(.0423xfuel burned) in slug-ft 


ynamic Coefficient Equations: 

Cd=Cdo +.5 x Cl x sin( &) 

Cdo= -.1025 x MACH + .5015 (1.8< MACH < 2.2) 
Cm=Cl (Xeg-Xcp)/d 


C1=Cl x og 


Key Wind Tunnel Results (per degree): 


MACH=1.8 


Cdo=.317 Cmq=-203 


« fe] 2 4 6 8 10 
cl ce] .18 .38 -62 -92 1.3 
Xcp 45.4 45.6 48.2 49.7 S0.2 50.0 
MACH=2.2 
Cdo=.276 Cmq=-190 
o fe) 2 4 6 8 10 
cl 0 19 .40 .65 .95 1.3 
Xcp 42.9 42.9 44.6 45.1 46.0 44.6 
Cmq=d¢(Cm) /d(qd/2v) 
Initial Muzzle Spin Rate 30 PI radians /sec 


in 


dee 


APPENDIX B TRANSIENT ANGLE OF ATTACK CALCULATIONS 


The following equation of motion applies to the dynamics 
of angle of attack (alpha) oscillations due to the torque 


equation deacribed in the three-degrees-of-freedom model. 


ny = Lw which is proportional to alphact) 


where angle of attack is of the form ; 


-§t 


e<(t)=e sintwt) where Sis the damping constant 


and w is the angular velocity. Therefore: 


tT 
5 x(t) dt dT 
oOo 


UT 
5 N e-St sin(wt) dt dT 
(dk =) 


we 


my 


dT 


= wt(§2 +w2) -28w + 
e“St(S2siniwt) +28 weostwt) -w2sin(wt) 


C82 + w2)2 


When t=0, the equation reduces to: 

-2§w +1(28w) = O which is expected if the initial 
condition is zero.(ie. dispersion at zero equals zero at 
time zero). 

As t approaches infinity the equation value approaches 
infinite dispersion, but the oscillations of alpha have 
damped out. Note that the first derivative, proportional to 


velocity, tends to a constant value: 


Cw/ (82 + w2 ))/m as t approaches infinity which is 


also to be expected. 


APPENDIX C-PSEUDOQ-VACUUM TRAJECTORY CONCEPT 


The basic concepts behind the pseudo-vacuum are not 
always readily apparent.This appendix analyzes the basic 
equationa of motion in the following sequence of cases to 
dlluatratea the concept: 

vacuum trajectory 
atmosphere trajectory 
thrusted trajectory 
wind influenced trajectory 
thrust equal drag trajectory 
In Figure C-1 (curve 1) the applicable equations of motion 


are: - 
mX1 ce) 


ee 
M21 “ag 


An accelometer on the projectile in the 2 direction would 
measure no acceleration as it is in free fall already. This 
is similar to a satellite orbit or someone on an elevator in 
free fall. 


When atmosphere is added the key equations are: 


m¥2 = -Dy = -1/2 9 Uy2 A Ca 


mZ2 -Dz - mg 


where U is relative wind and D is atmospheric drag.In this 


case Ux=Vx and Uz=Vz.This trajectory ia represented by curve 


2.An accelometer in the x direction measures Dx/m and Dz/n 
in the z direction. 
If thrust is added, the trajectory is as shown at curve 


3(thrust greater than drag).The key equationa are then: 


mx3 = -Dy + Ty 
MZ = -De + Tz “mg : 
' 
Accelerometers would measure (Tx-Dx)/m and (Ty-Dy) /m 


respectively.If thrust equals drag the trajectory would be 


curve 1. 


Vw 


Figure C-1 Pseudo-Vacuum Trajectory Concept Curves 


893 


In Figure C-1 wind (Vw) has been added for the curve 
in the x direction.The new equations are : 


mxq = (-Dy + Dw) + Tx 


mz4 -Dz + Tz “Mg 


where Dx + Dw= 1/2 9 Uy? Cg A 


and Ux(wind relative to the projectile) =Vx+Vw. 
thrust curve 4 results. 
Finally, if thrust is added to 
following equations are valid: 
mx5 = O or Tx = Dy + Dw 


RZS -mg or Tz = Dz 


The accelerometer(s) read the acceleration caused by any 
differences in the thrust and drag forces (zero in this 
case). 

It is clear that curves 1 and 5 are equal (and 3 if 
thrust=drag) but the burntime of a ramjet will not be equal 


in each case due to different drag profiles. 


t = (masa fuel) (Isp)/thruat 


Therefore headwinds will require more thrust and 


burntime will be shorter.Conversely tailwinds allow longer 


burning.Burnout will occur on the vacuum trajectory but 


at the same point, as shown in Figure C-2. 


not 


Figure C-2 Burnout Variations 
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APPENDIX D PSEUDO-VACUUM TRAJECTORY PROGRAM 


18 DIM X(300),Y (380), RO(ZQQ), D( SAN), MDOT (SAM), M( Baa, visa) 

2 DIM CDO(3aa), MA‘ 3aa) 

1S DIM T(3Zaa) 

20 MMIN=198 

30 ISPp=95a 

4@ DELT=.5 

SQ M(1) #225 

62 THETA=. 78539816340020068 

7@ VO=2e452 

75 LPRINT "TIME", "xX", "Y", "MASS", "SPEED" 

Q8@ FOR I=1 TO saa 

90 X(I)=VO*COS (THETA) *T (I) 

100 Y(1)=(VO*SIN(THETA) *T (I) )- (16. Q87*T (1) 2) 

110 V(1L)=( (VO*COS (THETA) ) *2+ (VO*SIN (THETA) — 
32.174*"T(1))*2)%.5 

120 TEMP =518. 688-. 2aS5662*Y (1) 

15@ RO(I)=(2116. 22% (1-(. QVZSEG2*Y (1) /518. 688) )*5. 256) / — 
(1716. 5*TEMP) 

155 IF Y(1I))=36089! THEN RO(I)=.@01271 

157 IF Y(I))=36089! THEN TEMP=389. 99 

158 MA(I)=V(I)/SGR(1. 4#1716. 3* TEMP) 

159 CDO(1)=.5@15-. 1025*MA(I) 

162 D(I)=. 276%. 5#RO(1) *. 3491*V(1) 42 ; 

161 D(I)=D(1)/. 276*CDO(I) 

170 MDOT(I+1)=D(1I)/ISP 

189) M(I+1)=M(1)—MDOT (1) *DELT 

190 T(I+1)=T(1)+DELT 

195 LPRINT TCI), X(I),Y(I),M(I),VCI) 

200 IF M(I+1) (=MMIN GOTO z2a 

219 NEXT I 

220 LPRINT "RAMJET BURNOUT BETWEEN T="sT(I) s"AND"sT(I+1) 
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REPROOUCED AT GOVERNMENT EXPE’ 


TIME x Y MASS SPEED 
a Q a 225 2452 
5 866.913 862.891 225 2442. 65 
1 1733. a3 1717.74 224. 637 2429. 36 
1.5 2600. 74 2564. 54 224. 286 2418.11 2 
2 3467.65 2403.3 223.947 2406. 93 ry 
2.5 4334.57 4234. 2 223.618 2395.8 
3 5201.48 5256.7 223.3 2384.73 
3.5 6268. 39 5871.32 222.992 2373.71 
4 6935.3 6677.91 222.695 2362.75 
4.5 7a@2. 22 7476. 46 222. 487 2351.85 
5 8669.13 8266. 36 222.128 2341.01 
5.5 9536. 04 9049.41 221.858 2250. 24 
6 10403 9823. 82 221.597 2319.52 > 
6.5 11269.9 1059.2 221.344 2528. 86 

_ 7 12136.8 11348.5 221.2399 2298. 27 
7.5 13203.7 12098.8 220. 861 2287.74 
) 1387@.6 12841 220.631 2277.28 
8.5 14737.5 13575.2 220. 429 2266. a9 
=) 15604.4 14301.4 228.193 2256. 55 
9.5 16471.3 15019.5 219.983 2246.29 
12 17338. 2 15729.6 219.781 2226.1 » 
13.5 18205. 2 16431.6 219. 584 2225. 98 
11 19872. 1 17125.6 219.394 2215. 92 
11.5 19939 17811.5 219. 209 2205. 94 
12 20805. 9 18489. 4 219.03 2196. 33 
12.5 21672.8 19159.2 218. 856 2186.2 
13 22539.7 19821 218.687 2176. 43 . 
13.5 23426. 7 20474.8 218.523 2166.75 Sen 
14 24273.6 21120.5 ° 218. 364 2157.14 » 
14.5 25148.5 21758.2 218.21 2147.61 
1s 26207. 4 22387.8 218. 261 2138.15 
15.5 26874.3 23009. 4 217.915 2128.78 
16 27741.2 23622. 9 217.774 2119.49 
16.5 28608. 1 24228.4 217.637 2118.2 5 
17 29475 24825.9 217.504 2101.14 Shee 
17.5 20342 25415.3 217.374 2092.1 ideas 
18 31208.9 25396.7 217.248 2083. 14 , 
18.5 32075.38 26570 217.126 2074.27 
19 32942. 7 27135.3 217. 087 ZOES. 48 
19.5 3289.6 27692.5 216.891 2056. 78 
2a 34676.5 28241.7 216.778 2048. 17 
20.5 35543.4 28782.9 216. 669 2039.65 
21 36419.3 29316 216. S62 2331.23 
21.5 37277.3 29841 216. 458 Sde2. 89 
22 38144.2 30258.1 216. 357 2014. 65 ’ 
22.5 3911.1 30867 216.258 2026. 51 
22 39878 313638 216. 162 1998. 46 
22.5 40744, 9 31860. 9 216. 268 1992. 51 
24 41611.8 32345.7 215.977 1982. 66 
24.5 42478.7 328e2.5 215. 888 1974.9 
25 43345.7 33291.3 215. 8@1 1967.25 
25.5 44212. 6 33752 215.717 1959.7 
26 4579.5 34204. 7 215.634 1952. 26 » 
26.5 45946. 4 34649.3 215.553 1944.91 
27 46812.3 25085. 9 215.474 1937. 68 
27.5 47688. 2 35514.4 215.397 1932.55 
2a 48547. 1 25934.9 215. 222 1923.53 
23.5 49414 36347.4 215.249 1916.62 
29 Sv2a1 36751.8 215.177 1909. Be 
23.5 S1147.93 S7148.¢ 215.848 L9BsZ1F 
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